Spontaneous spirals in vibrated granular chains 
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We present experimental measurements on the spontaneous formation of compact spiral struc- 
tures in vertically-vibrated granular chains. Under weak vibration, when the chain is quasi two- 
dimensional and self- avoiding, spiral structures emerge from generic initial configurations. We com- 
pare geometrical characteristics of the spiral with that of an ideal tight spiral. Globally, the spiral 
undergoes a slow rotation such that to keep itself wound, while internally, fast vibrational modes 
are excited along the backbone with transverse oscillations dominating over longitudinal ones. Spi- 
rals have an extremely small volume in phase space, and hence, their formation demonstrates how 
nonequilibrium dynamics can result in a nonuniform sampling of phase space. 

PACS numbers: 81.05. Rm,83.10Nn 



I. INTRODUCTION 

Ratchet mechanisms are thought to play an important 
role in producing net transport or motion by asymmet- 
rically absorbing energy from a stochastic source [1-4]. 
Such mechanisms should apply to systems far from equi- 
librium as well, specifically with respect to the formation 
of patterns or structures. In other words, spontaneously 
generated patterns might arise from a stochastic forcing. 
Shaken granular media are natural systems for investi- 
gating possible ratchet effects in nonequilibrium systems 
[5]. Mechanical vibration is athermal and coupled with 
the dissipative collisions results in effectively a stochastic 
forcing often manifested by anomalous velocity fluctua- 
tions [6-8]. 

Recently, we have explored a chain of coupled hollow 
spheres attached by rods and vibrated vertically on a 
horizontal surface [9,10]. The vertical vibration pumps 
energy into the chain and, because the collisions with the 
plate are typically asynchronous with the drive, the chain 
develops irregular vertical structure. Collisions between 
balls on the chain generate a noisy dynamics of lateral 
motions so that viewed from above the chain appears as 
effectively two-dimensional. Further, when the vibration 
is weak so that the vertical displacement of the chain is 
less than one ball diameter, the chain cannot cross itself 
and is self-avoiding. 

Previous experiments suggest that vibration induces 
strong fluctuations which in turn cause the chain to ex- 
plore all regions of the available phase space, similar to 
the ergodicity principle in equilibrium statistical mechan- 
ics. Statistical properties of the unraveling times of knots 
is consistent with random thermal motion of the cross- 
ing points constituting the knot [9]. Moreover, two-loop 
figure-eight configurations in a ring chain prefer states 
with maximal disparity in the size of the loops, consis- 
tent with uniform sampling of the microscopic configura- 
tion space [10]. Further, the statistical properties of free 
and closed chains [11,12] are reminiscent of the proper- 
ties of ideal polymers in equilibrium systems [13,14] for 
such quantities as, for example, the radius of gyration. 



One would expect on the basis of these previous results 
that a free chain weakly excited so that it cannot cross it- 
self would explore phase space uniformly. Such is not the 
case. Instead the chain may form a compact spiral, see 
Figure la, starting from some generic initial condition. 
From the available phase space, the system condenses 
into a very small region, namely the spiral state. In this 
paper, we present a characterization of the spiral state. 
We start with a short description of the experimental 
apparatus, followed by discussions of the formation, ge- 
ometry and dynamics of the spirals. 
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FIG. 1. The raw image (a) and the reconstructed spiral (b). 
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II. EXPERIMENTAL SET-UP 

The experiment consists of a vertically-vibrated 
anodized-aluminum plate driven harmonically by an elec- 
tromagnetic shaker. The plate was circular with a diame- 
ter of 11.50 inches. It was coupled to the shaker by a shaft 

that occupied the bore of a precise square-section linear 
air bearing. As a result, the plate's vertical vibration 
was accurately sinusoidal with amplitude A and angular 
velocity lo = 27r/ where / is the driving frequency. In ad- 
dition, the square-section of the bearing ensured that the 
plate did not experience rotational torques. An acrylic 
lateral boundary confined the chain to the plate. The 
instantaneous acceleration of the plate, measured with 
an accelerometer, was maintained at constant value of 
the dimensionless peak acceleration T = a/g = (jJ^A/g to 
within 2%. 

The commercially-available ball chain consisted of N 
hollow, approximately spherical shells that were con- 
nected by dumb-bell shaped rods. The chains were 
nickel-plated steel with a nominal bead diameter d = 
inches with a variation of about 1%. The rods had a di- 
ameter 0.22(i and constrained the inter-bead separation 
1 < & < 1.4 measiued in units of the bead diameter. The 
maximal angle between two consecutive rods was w 50° 
so that the smallest circle that could be made from the 
chain had n = 8 beads. The ratio of the diameter of the 
available plate to the bead diameter is about 110, and 
boundary effects were negligible. 

In our experiments, a ball chain was placed onto the 
vibrating plate. The relevant driving conditions are the 
frequency, /, and the acceleration, T . We restricted our 
attention to weak vibration, F < 2, where the geome- 
try is quasi two-dimensional. The chain is initiated in a 
random configuration subject only to the condition that 
it did not intersect itself. In these weak driving condi- 
tions, the chain remains close to the plate surface, free of 
self-intersections . 

We obtained digital images of the chain using a high- 
resolution (1024 x 1024 pixels) CCD camera. Data sets 
were acquired at frame rates of 1 Hz to observe the slow 
global dynamics and at 250 Hz to characterize internal 
vibrational modes. Using standard image analysis proce- 
dures, we extracted positions of the bead centers. These 
were then ordered to recover the locus of the spiral. For 
more details see [10]. A typical image and chain recon- 
struction are shown in Figure 1. 



III. SPIRAL FORMATION 

We find experimentally that sufficiently long chains 
spontaneously organize into spirals in a narrow range 
of accelerations, 1.70 < F < 1.85, and frequencies 
10 < / < 25 Hz. Starting from a linear chain, a tight 
spiral core nucleates (either clockwise or counter clock- 
wise) at one end of the chain. This core rotates, and 



as a result the chain winds into a spiral. The formation 
of compact spirals from a nucleated core typically takes 
3—7 minutes for a chain with N — 256 beads. However, 
the nucleation of the spiral core from a linear chain is a 
random event occurring on the scale of 30 minutes. Once 
formed, the spiral was stable. 

The formation of spirals is more sensitive to variations 
in the acceleration rather than in the frequency, similar 
to pattern formation and clustering in vibrated granu- 
lar layers [15,16]. When the acceleration is too large, the 
chain exhibits strong excitations that quickly destroy the 
spiral. On the other hand, when the acceleration is too 
small, lateral chain motion is insufficient to nucleate the 
spiral core. 

Interestingly, for a given chain, the spiral always nucle- 
ates with its core at the same end. Moreover, artificially 
formed spirals with their core at the opposite end are 
unstable. The mechanism responsible for this symme- 
try breaking are minute imperfections in size, shape, and 
weight of the extremal beads. Asymmetry underlies nu- 
merous transport phenomena in physical and biological 
systems, for example, Brownian motors [3]. Typically, 
asymmetry in the environment coupled with oscillatory 
driving leads to ratchet motion. While similar in spirit, 
granular chains are different in that the asymmetry is in 
the transported object itself. Such a transport mecha- 
nism may be relevant in biological systems, for example, 
it can be realized via defects in the macromolecules. 



IV. SPIRAL GEOMETRY 

The asymmetries are important in the nucleation of the 
spiral core and its evolution to a compact rotating state. 
Since we are unable to systematically control these asym- 
metries, we do not consider here the transient formation 
stage, but rather, we focus on the characteristics of the 
formed spiral. As a reference, we use a ideal tight spiral. 
Consider such a spiral formed by a flexible rope. With- 
out loss of generality, its diameter is set to unity. A point 
along the spiral is (over)specified by its arc length s, ra- 
dius r, and winding angle 0. It is convenient to think of 
the spiral as composed of concentric tori. In such an ar- 
rangement, the addition of the outermost torus increases 
the angle by /S.6 = 27r and the arc length by As = 27rr. 
Since Ar = 1 we have ds/dr = 27rr and ds/d6 = r. It 
follows that r(s) = s/sJ-k and 9{s) = v^Itts. 

Consider now a bead chain. The bead number n is a 
natural measure of the arc length, and assuming an aver- 
age separation of (6), one has s(n) = n{b). Moreover, the 
distance between the nth bead and the 1st bead, r(n), is 
used as a measure of the radius, since the spiral center is 
not uniquely defined. The (s, 9) coordinates of the first 
bead are (1, 0). Substituting s(n) = n{h) in the above ex- 
pressions for 9{s) and r(s), and using the law of cosines 
yields an expression for the winding angle and distance 
to the core as a function of the bead number: 
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e{n) = V47rn(6), 



(1) 



r(n) 



l + !^-2jr^cose{n). (2) 



Of course, these expressions apply to sufficiently large 
spirals and 1. 

We have measured the average angle and distance as 
a function of the bead number from over 500 frames. In 
these experiments and for the remainder of this paper we 
have N = 256 beads, / = 16 Hz and T = 1.77. In re- 
constructing the spiral from the image our ordering was 
such that the innermost core bead corresponds to n = 1. 




A 

B 

V 



60 



50 



40 



30 



20 



"T 1 1 1 1 r 



(c) 



10 



_l I I I I L_ 



1 



2 3 4 5 6 7 

k 



FIG. 2. Tiic mean angle {6)/2n (a) and distance (r) (b) as 
a function of ^/n. In (c) the mean number of beads (m) per 
turn is plotted vs the turn number k. 

In Figures 2a and b we have plotted the averages 
{6){n)/2'K and {r){n) versus ^/n respectively. For com- 



parison, wc also show the tight spiral angle and distance 
calculated from Equations (1) and (2). Hero, the aver- 
age separation (6) = 1.2 was measured directly, as dis- 
cussed below. The slower (faster) linear growth of {9) (n) 
((r) (n)) corresponds to the spiral becoming less tight the 
further one goes from the compact core. Nevertheless, 
the gross features such as the overall linear growth of 
(r(n)) with ,Jn, and the periodic corrections, with a pe- 
riod that is fixed in ^/n, indicate that the spiral is ap- 
proximately tight. 

The periodic oscillations of {r{n)) allows for a char- 
acterization of the number of beads in a given turn of 
the spiral. A single turn is contained between two con- 
secutive minima in Figure 2b. As shown in Figure 2c, 
the mean number of beads per turn grows linearly with 
the loop index, consistent with the expected behavior for 
concentric tori. The deviation from linearity for the last 
turn indicates that the spiral has a straight tail. From 
this measurement, the average size of the tail is estimated 
at 10, consistent with direct observation. We also find 
that different sets of 500 images, corresponding to inde- 
pendently formed spirals, are characterized by quantita- 
tively similar mean angle, distance and number of beads 
per turn as in Figure 2. 

When comparing with a tight spiral, we implicitly as- 
sumed that the separation between two beads is constant. 
Direct measurement shows that this is a reasonable as- 
sumption as variations in (6) are of the order of 10%. 
However, this deviation has a clear trend. The separa- 
tion increases as a function of the bead number and it 
saturates at a value of (6) w 1.24 for n > 150. 
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FIG. 3. The average bead separation (b) (measured in units 
of bead diameter d) as a function of the bead number n. The 
thick line is a 12-bead running average. 



V. SPIRAL DYNAMICS 

Once formed, the vibrating spiral has dynamics on fast 
and slow time scales. In Figure 4a, we have plotted an 
overlay of 16 consecutive reconstructions of the spiral ob- 
tained at a sampling rate of 250 Hz. The 16 images 
correspond to approximately 1 oscillation of the plate. 
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The lateral wave modes visible here constitute in part 
the fast dynamics of the spiral. What is not obvious in 
Figure 4a is the large number of collisions and the con- 
strained motion that effectively gives rise to the trans- 
verse modes. Approximately 4 beads separate the nodes 
of these modes. On a much longer time scale we find that 
the spiral undergoes axial rotation. 




FIG. 4. Sixteen consecutive reconstructions (a) of a spiral 

during 1 oscillation showing the fast transverse wave modes. 
Images with a separation of 15 seconds (b) illustrating the 
slow rotation that keeps the spiral wound. 

Four snapshots of a spiral with a separation of 15 sec- 
onds or 240 oscillations between snapshots is shown in 
Figure 4b. The spiral rotates so as to keep tightly wound 
at a constant rate of about one revolution per minute. 
On shorter time scales the rotation is not constant, and 
is better described as a ratcheting motion with the spiral 
core moving back and forth. The asymmetry between 
the winding and unwinding motion of the core is trans- 
mitted along the spiral length and over long time scales 
produces the global rotation. 

We now focus on the richer fast dynamics. Denot- 
ing the displacement in position between consecutive im- 
ages for each bead S = {Sx,Sy), we find that they are 
independent of the bead position with slight differences 
for the end beads. We sampled in excess of 5000 dis- 



placements per bead. The data was obtained in sets of 
approximately 540 frame to frame displacements per ex- 
periment sampled at 250 Hz. We verified that for each 
bead the mean value (S) vanished and that the displace- 
ments were isotropic. Since the spiral rotates, its fast dis- 
placements sample all directions in the plane so that the 
orientation of the coordinate axes is irrelevant. We have 
plotted the probability distribution functions P{Ax) and 
P{Ay) in Figure 5. Here, we have defined the normal- 
ized displacements A^^y = Sx^y/cr^^y where a^^y are the 
respective standard deviations. Our results show that the 
quantities and are identically distributed with sig- 
nificant incidence for displacements as large as ±3ax,y 
The probability distribution functions are symmetric and 
non-Gaussian. Over populated (with respect to a Gaus- 
sian) velocity distributions [7] and, in particular, expo- 
nential tails [8,17] are actually observed in excited gran- 
ular gases and are understood to be a consequence of 
the dissipative bead-bead and bead-plate collisions. Our 
observations suggest that velocity statistics of beads in 
granular chains may be anomalous as well. Although the 
exponential tail of the displacement distribution suggests 
that the velocity distribution has a stretched exponential 
tail, we are unable to discern its precise form due to in- 
sufficient time resolution. 
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FIG. 5. The probability distribution functions of the nor- 
malized displacements A^, and Ay for all beads of the spiral. 

At present we only observe the lateral fast dynamics. 
The transverse wave modes illustrated in Figure 4b also 
appear in the vertical plane. Although the in-plane mo- 
tion consists primarily of transverse modes, they must 
be accompanied by longitudinal modes that allow for the 
stretching and contraction of tlw^ spiral. To further study 
these modes we consider the locus of the spiral given by 
the bead positions Xj = {xi,yi) and inter-connecting rods 
r, = Xj+i— Xj. The mean locus of the spiral (r)j is defined 
as the time average of (t) over an experimental run. For 
the frame acquired at time t, we calculate the spiral dis- 
placements Ar i{t) = ri{t) — (r)j. The longitudinal Li{t) 
and transverse Ti{t) displacements from the mean spiral 
arc defined by Ar,{t) • (r),/|(r),| and Ar i{t) x (r)j/|(r),|, 
respectively. In Figure 6a we have plotted Ti{t) as a func- 
tion of position and time. The space-time plot shows that 
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between collisions with the plate, positive (negative) dis- 
placements continue as positive (negative) i.e. the beads 
are roughly ballistic between collisions. In Figure 6b 
we plot representative displacements for 80 oscillations 
of the plate. We see that the transverse displacements 
are roughly periodic and are significantly larger than the 
longitudinal ones. This behavior is typical for all inter- 
connecting rods. 
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FIG. 6. A space-time plot (a) of the transverse displace- 
ments of the spiral. Position (horizontal axis) along the spiral 
backbone is ordered from core to tail and time (vertical axis) 

spans 64 consecutive images corresponding to 4 drive cycles. 
Black (white) denote a positive (negative) displacement. In 
(b) Ti{t) (solid line) and Li{t) (dashed line) are plotted for 
i = 64. 



VI. CONCLUSION 

A general expectation from previous experimental 
work on granular chains [9-12] is that configurationally 
the chain assumes very many different states and in this 
manner explores a large region of its phase space. Con- 
trary to this, we found that a long chain would sponta- 
neously nucleate into a spiral in a narrow range of weak 
accelerations. It is surprising that a linear chain initiated 
in a random configuration invariably evolves to a highly 
specialized state that occupies a very small volume of the 
available phase space. Moreover, it is interesting that mi- 



nor differences between the two ends of the chain have 
such a pronounced effect. 

The origins of spiral formation in vibrated chains is the 
slight forward motion of the chain induced by variations 
at the ends of the chain. As a result of this asymme- 
try, the chain crawls forward until it encounters itself. 
It then has two choices: to turn inwards and wrap itself 
into a spiral or outwards forming locally aligned chain 
segments. In the latter case, the end eventually escapes. 
With successive encounters, it becomes likely that the 
end will be captured and coil up into a spiral. The rate of 
nucleation depends, therefore, on the forward ratcheting 
motion coupled with the possibility of becoming trapped. 
Such a mechanism would apply to any flexible forward- 
translating object confined to two- dimensions and mov- 
ing stochastically. 

Once formed, the spiral geometry and dynamics is very 
reproducible. The spiral becomes looser away from the 
core and eventually, its tail is practically linear. Dynam- 
ically, bead-bead and bead-plate collisions induce trans- 
verse and longitudinal modes, the former being stronger. 
On a longer time scale, the spiral rotates in the direction 
of coiling and so remains stable. 

It will be interesting to further study the spiral for- 
mation mechanism. This may be done for example by 
manipulating one of the chain ends using say a relatively 
larger bead. Observing the rotation rate as a function 
of the bead size may elucidate the physical nature of the 
ratcheting mechanism. Moreover, the nucleation process 
is in itself intriguing. The distribution of waiting times 
and its extremal characteristics may alternatively shed 
light on the same question. 

We thank Matthew Hastings, Charles Reichhardt, and 
James Glazier for useful discussions. This research is 
supported by the U.S. DOE (W-7405-ENG-36). 
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